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PURPOSE:  The  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  described  herein 
provides  guidance  on  coupling  the  regional  circulation  model,  Advanced  Circulation  (ADCIRC), 
with  the  local  circulation  model  M2D.  Coupling  of  models  is  an  efficient  and  accurate  means  of 
calculating  water  level  and  current  on  a  local  domain,  such  as  at  and  around  an  inlet.  This  CHETN 
also  serves  as  a  tutorial  for  coupling  the  two  models  within  the  Surface-Water  Modeling  System 
(SMS). 

BACKGROUND:  Circulation  modeling  at  tidal  inlets  requires  accurate  boundary  conditions  at 
the  seaward  edge  of  the  domain.  Regional  models  developed  by  the  Coastal  Inlets  Research 
Program  (CIRP)  provide  information  that  can  be  applied  as  boundary  conditions  for  local  (project 
scale)  models.  Examples  of  regional  models  available  from  the  CIRP  are  Long  Island,  NY, 
southwest  Washington,  northern  California,  east  coast  of  Florida,  and  the  Maryland  and  Delaware 
coast.  Water  levels  from  solutions  of  these  regional  models  can  be  applied  as  boundary  conditions 
for  local  projects. 

Large-scale  physical  oceanographic  processes  such  as  the  tide  and  storms  can  be  represented  in 
regional  models  and  the  information  cascaded  locally  by  furnishing  boundary  (driving)  conditions  to 
meet  project  analysis  requirements.  Also,  multiple  projects  with  great  detail  can  be  included  within 
a  computationally  efficient  regional  domain  (Militello,  Kraus,  and  Brown  2000).  Benefits  include 
automatic  calculation  of  mutual  interactions  among  projects,  such  as  multiple  inlets  sharing  the  same 
bay  or  lagoonal  system,  and  the  availability  of  baseline  information  from  which  new  projects  can  be 
added.  In  addition,  duplication  of  effort  can  be  avoided,  increasing  efficiency  and  cost-effectiveness. 

Engineering  applications  require  numerous  simulations  with  modified  grids  to  determine  the 
consequences  of  project  alternatives.  An  efficient  means  of  conducting  multiple  simulations  is  to 
apply  a  local  model,  which  eliminates  the  computational  overhead  of  a  large  domain.  By  coupling 
regional  and  local  models,  the  local  model  can  propagate  hydrodynamic  conditions  on  the 
boundaries  calculated  by  the  regional  model  into  its  smaller  domain.  Coupling  of  the  models  is 
accomplished  by  passing  time  series  of  water  level  obtained  from  the  regional  model  solution  to  the 
local  model  as  boundary  conditions.  Because  the  boundary  conditions  for  the  local  model  were 
calculated  on  a  regional  domain,  they  contain  spatial  variations  that  preserve  such  features  as  tidal 
phase,  surge  distribution,  and  wind  setup  or  setdown.  Other  boundary  conditions  or  input,  such  as 
river  discharge  or  wind,  can  be  specified  for  the  local  model  as  required  by  the  project. 

The  coupled  modeling  approach  is  advantageous  because  the  local  model  can  apply  ADCIRC- 
calculated  water  levels  from  previous  studies  as  boundary  conditions  if  they  are  available  for  the 
time  period  of  interest.  Alternatively,  ADCIRC  can  be  run  for  the  time  intervals  required  for  the 
study,  and  then  the  local  model  can  be  applied.  In  both  situations,  modifications  to  the  local  grid  can 
be  easily  made  for  project  alternatives  without  the  necessity  of  multiple  simulations  with  the  large 
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regional  model.  This  approach  minimizes  the  time  required  for  the  user  to  work  with  the  regional 
finite-element  model. 

This  CHETN  describes  coupling  of  the  regional  oceanographic  model  ADCIRC  (Luettich, 
Westerink,  and  Scheffner  1992)  with  the  local  model  M2D  (Militello,  Zundel,  and  Buonaiuto,  in 
preparation;  Militello  1998).  Automated  coupling  of  the  models  is  implemented  within  the  SMS 
(Zundel  2002)  through  interpolation  of  ADCIRC-calculated  water  levels  to  M2D  boundary  cells 
(Militello  and  Zundel,  in  preparation).  In  the  following,  the  circulation  models  are  introduced  and 
information  on  coupling  them  within  the  SMS  is  provided.  An  example  of  automated  coupling  is 
then  given. 

ADCIRC:  The  two-dimensional  depth-integrated,  finite-element  circulation  model  ADCIRC  was 
developed  for  large-scale  applications  within  which  fine  detail  can  be  specified  for  local  areas. 
ADCIRC  has  great  flexibility  in  its  element  sizes,  which  is  convenient  for  creating  calculation 
meshes  in  which  coarser  resolution  is  specified  in  areas  away  from  the  study  site  to  improve 
computation  time.  Resolution  can  be  increased  in  areas  of  strong  bathymetric  gradients,  such  as  over 
the  continental  shelf,  to  maintain  computational  accuracy. 

Features  of  ADCIRC  include  multiple  forcing  options,  robust  wetting  and  drying,  hot  start 
capability,  calculation  in  spherical  or  Cartesian  coordinates,  choice  of  bottom-stress  models,  internal 
harmonic  analysis,  and  user-specified  global  and  station  output.  Forcing  can  be  imposed  by  tidal 
constituents,  water-surface  elevation  time  series,  flow-rate  time  series,  wind  stress,  atmospheric 
pressure,  and  wave  stress.  Information  on  these  and  other  model  features  are  provided  in  an  online 
ADCIRC  User’s  Guide  that  can  be  accessed  on  the  internet  at  http://www.unc.edu/depts/ 
marine/C _C ATS/adcirc/document/ ADCIRC _main  Jrame.html.  The  Web  site  also  lists  ADCIRC- 
related  publications  on  research  and  model  applications. 

Through  the  large-domain  capability  of  ADCIRC,  processes  on  multiple  spatial  and  temporal  scales 
can  be  calculated.  For  example,  nonlinear  tidal  propagation  over  the  continental  shelf  is  a  large- 
scale  process  that  can  be  combined  in  the  same  simulation  as  flow  though  a  tidal  inlet.  Calculated 
water  level  and  current  velocity  values  represent  the  response  of  the  combined  processes  on  all 
scales  included  in  the  simulation.  Thus,  transfer  of  water  levels  calculated  with  ADCIRC  to  a  local- 
domain  model  preserves  the  description  of  the  processes  at  the  boundaries  of  the  local  model. 

M2D:  The  two-dimensional  depth-integrated  finite-difference  hydrodynamic  model  M2D  was 

developed  for  project-scale  applications,  such  as  inlets,  navigation  channels,  and  bays.  M2D 
operates  on  an  orthogonal  rectilinear  grid  in  which  cell  side  lengths  can  vary.  Greater  resolution  can 
be  specified  in  areas  where  it  is  needed,  such  as  at  structures  and  areas  of  changing  depth  (channels, 
shoals,  wetlands  and  marshes),  and  special  target  areas  of  interest.  In  addition,  the  grid  can  be 
aligned  as  needed  for  a  specific  application.  For  example,  for  an  inlet  application,  the  grid  can  be 
aligned  with  its  x  or  y  axis  oriented  along  the  principal  axis  of  the  channel.  This  capability  for  grid 
rotation  provides  flexibility  for  setting  the  orientation  to  maximize  efficiency  and  accuracy. 

Features  of  M2D  include  robust  flooding  and  drying,  wave-stress  forcing,  wind-speed  dependent 
(time-varying)  wind-drag  coefficient,  variably-spaced  bottom-friction  coefficient,  efficient  grid 
storage  in  memory,  hot  start  capability,  and  user-specified  global  and  station  output.  M2D  can  be 
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forced  with  any  combination  of  water  level,  tidal  constituents,  flow  rate,  wave  stresses,  and  wind. 
Model  output  is  written  in  formats  convenient  for  importing  into  commercially  available  plotting 
packages.  The  rectilinear  grid  of  M2D  is  convenient  to  create  and  modify  as  compared  to  a  finite- 
element  (irregular)  mesh,  especially  for  those  with  modest  experience  in  modeling  or  for  fast-track 
projects.  Also,  sediment- transport  and  morphology  change  algorithms,  presently  under  development 
in  the  CIRP,  are  readily  and  naturally  accommodated  on  a  rectilinear  grid. 

MODEL  COUPLING:  Coupling  of  ADCIRC  and  M2D  takes  place  through  water-level  forcing 
boundary  conditions  specified  for  M2D.  Time  series  of  water  level  calculated  by  ADCIRC  serve  as 
forcing  information  for  M2D.  The  basic  procedure  for  model  coupling  is: 

a.  Run  ADCIRC  model  for  simulation  interval. 

b.  Specify  ADCIRC  water-level  files  as  boundary  forcing  for  M2D. 

c.  Extract  water-level  values  from  ADCIRC  solution  to  be  applied  as  forcing  for  M2D. 

d.  Run  M2D  for  simulation  interval. 

Coupling  of  the  two  models  is  automated  within  the  SMS.  Specification  of  ADCIRC  water-level 
files  and  extraction  of  values  for  M2D  are  conducted  based  on  an  existing  ADCIRC  simulation  and 
the  boundary  condition  settings  for  M2D.  The  SMS  interpolates  the  ADCIRC-calculated  values  to 
M2D  cell  centers  where  water  level  is  specified. 

For  accurate  calculation  of  boundary  values  for  M2D,  the  ADCIRC  mesh  must  contain  adequate 
resolution  near  the  M2D  forcing  boundary.  If  the  ADCIRC  mesh  resolution  is  too  coarse,  then  the 
interpolation  of  water-level  values  to  the  M2D  boundary  cells  may  be  inaccurate  or  biased, 
depending  on  the  placement  of  the  ADCIRC  computational  nodes  relative  to  the  M2D  cell  centers. 
Appropriate  resolution  of  the  ADCIRC  mesh  will  depend  on  the  domain  and  requirements  of  a 
study.  Thus,  the  user  must  evaluate  what  resolution  is  necessary. 

Figure  1  shows  an  ADCIRC  mesh  at  an  M2D  boundary.  The  models  were  coupled  for  a  hydro- 
dynamic  and  sediment  transport  study  at  Bay  Center,  WA  (Militello  2002).  ADCIRC  computational 
nodes  are  located  at  the  vertices  of  the  triangles.  Water-level  forcing  values  for  M2D  are  supplied  at 
the  centers  of  the  M2D  boundary  cells,  denoted  by  the  red  dots.  The  fine  resolution  of  the  ADCIRC 
mesh  relative  to  the  M2D  grid  in  this  example  was  deemed  necessary  because  there  was  significant 
variation  in  tidal  phase  across  the  boundary.  This  level  of  detail  in  the  ADCIRC  calculations 
preserved  the  tidal  phase  over  the  M2D  boundary,  resulting  in  accurate  calculations  on  the  local 
(M2D)  domain. 

AUTOMATED  COUPLING  WITHIN  SMS:  The  SMS  contains  model-specific  interfaces  for 
both  ADCIRC  and  M2D.  These  interfaces  provide  a  suite  of  tools  for  finite-difference  grid  and 
finite-element  mesh  development,  control  file  setup,  boundary  condition  specification,  and  grid/mesh 
quality  checks.  Runs  can  be  launched  from  within  the  SMS  and  results  loaded  for  post-processing 
and  visualization. 

To  take  advantage  of  the  efficiency  of  running  local  models,  SMS  includes  tools  to  simplify 
extraction  of  boundary  condition  data  from  the  regional  solution  at  the  local  model  boundary 
locations.  The  extraction  process  for  creation  of  boundary  conditions  for  the  local  model  relies  on 
the  existence  of  a  regional  model  with  a  solution.  Once  water-level  values  are  obtained  from  the 
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Figure  1 .  ADCIRC  nodes  and  M2D  cells  along  a  boundary 

regional  model,  SMS  interpolates  them  to  the  centers  of  the  M2D  forcing  boundary  cells.  In  setting 
up  a  simulation  for  the  local  model,  the  user  specifies  the  location  for  boundary  conditions  by  means 
of  a  cell  string.  A  cell  string  is  a  list  of  cells  where,  in  this  case,  the  water-surface  elevation  (WSE) 
will  be  forced.  The  boundary  condition  is  applied  to  a  cell  string  by  selecting  from  a  list  of 
supported  boundary  conditions,  as  shown  in  Figure  2. 


Figure  2.  M2D  Boundary  Conditions  dialog 


4 


ERDC/CHL  CHETN-IV-42 
June  2002 


The  SMS  detects  whether  a  regional  model  exists  in  memory,  and  it  then  activates  the  extract 
capability  as  a  possible  source  of  the  boundary  condition.  Then,  to  extract  a  boundary  condition,  the 
user  selects  the  Extract  BC  option.  This  action  displays  the  Extract  Boundary  Condition  Dialog, 
shown  in  Figure  3.  In  this  dialog,  the  user  selects  the  data  set  from  the  regional  model  supplying  the 
boundary  condition  in  the  local  model.  The  dialog  also  allows  the  user  to  specify  the  time  range  to 
be  extracted  from  the  regional  model  and  a  potential  offset  for  the  time  into  the  local  model. 
Conversion  of  time  units  is  also  administered  automatically.  It  should  be  noted  that  the  regional  and 
local  models  must  be  in  the  same  coordinate  system  for  the  extraction  and  interpolation  to  function 
correctly.  SMS  includes  tools  to  convert  information  from  one  coordinate  system  to  another.  If  the 
model  grids  have  been  computed  in  different  coordinate  systems,  the  conversion  can  be  done  by 
reading  in  the  regional  model  grid  and  converting  it  to  the  coordinate  system  of  the  local  model  grid 
before  the  local  model  is  brought  into  SMS. 


Figure  3.  Extract  Boundary  Condition  dialog 

Once  the  water-surface  elevation  curves  are  extracted,  they  can  be  viewed  individually  or  as  a  group 
with  the  View  Extract  button  in  the  M2D  Boundary  Conditions  dialog  (Figure  4).  All  water-surface 
elevation  curves  are  then  stored  with  the  M2D  input  data  as  part  of  the  local  model  run 
specifications. 

EXAMPLE  APPLICATION :  To  demonstrate  the  coupling  capability,  an  example  application  at 
Shinnecock  Inlet,  Long  Island,  NY,  is  provided.  Shinnecock  Inlet  is  the  easternmost  inlet  on  the 
south  shore  of  Long  Island.  Problems  faced  at  the  inlet  include  chronic  erosion  on  the  downdrift 
beach,  scour  near  the  jetty,  and  strong  ebb  current.  A  regional  ADCIRC  model  (Figure  5)  was 
developed  for  a  hydrodynamic  and  morphologic  study  of  Shinnecock  Inlet  (Militello  and  Kraus 
2001). 
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Figure  4.  View  Extracted  Boundary  Condition  dialog 


Figure  5.  Regional  ADCIRC  mesh  for  Long  Island 
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Global  water-surface  elevations  from  the  Militello  and  Kraus  (2001)  study  were  saved  every  0.5  hr. 
An  M2D  grid  was  developed  for  Shinnecock  Inlet  and  Bay,  as  well  as  a  portion  of  the  nearshore 
area.  Figure  6  shows  the  M2D  grid  overlaid  on  the  local  area  of  the  ADCIRC  mesh.  Water  levels 
from  the  global  ADCIRC  solution  file  furnish  forcing  for  the  seaward  M2D  boundary. 


Figure  6.  M2D  grid  shown  with  ADCIRC  mesh  at  Shinnecock  Bay  and  Inlet 

M2D  and  ADCIRC  must  be  in  the  same  coordinate  system  for  the  coupling.  Because  M2D  requires 
a  Cartesian  coordinate  system,  it  is  preferable  to  convert  the  ADCIRC  grid  from  spherical  to 
Cartesian  coordinates  (if  it  has  not  already  been  saved  in  the  same  coordinate  system  as  M2D).  If  a 
conversion  is  required,  the  easiest  approach  is  to  load  the  ADCIRC  grid  first  and  convert  it  to  the 
same  coordinate  system  that  the  M2D  grid  lies  in.  Then,  the  M2D  project  can  be  loaded  into  the 
SMS.  Alternatively,  ADCIRC  can  be  saved  in  the  required  Cartesian  coordinate  system  and  loaded 
into  SMS  either  before  or  after  M2D.  For  this  example,  ADCIRC  and  M2D  have  both  been  saved  in 
New  York  State  Plane,  Long  Island,  with  units  of  meters. 

An  ADCIRC  global  solution  file  must  be  loaded  into  the  SMS  before  water-level  values  can  be 
specified  at  the  M2D  boundary.  The  ADCIRC  solution  can  be  loaded  into  the  data  browser  in  the 
mesh  module.  Once  the  ADCIRC  solution  has  been  loaded,  coupling  of  M2D  with  ADCIRC  can 
take  place. 

To  extract  ADCIRC  water  levels  from  the  global  solution  to  apply  as  M2D  boundary  conditions,  the 
user  must  switch  from  the  SMS  mesh  module  to  the  Cartesian  grid  module.  A  cell  string  must  be 
defined  along  the  seaward  and  lateral  boundaries  where  the  extracted  water-level  values  will  be 
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applied.  To  create  the  cell  string,  click  on  the  Create  Cellstring  tool.  Next,  click  on  an  end  cell  of 
the  cell  string.  Holding  down  the  shift  key,  click  each  comer  along  the  boundary  until  all  boundary 
cells  are  included  in  the  cell  string,  and  press  Enter.  Figure  7  shows  the  newly-created  cell  string  in 
red. 


Figure  7.  M2D  grid  with  seaward  boundary  cell  string  shown  in  red 

The  next  step  will  define  the  cell  string  as  an  extracted  water-level  forcing  boundary.  Click  on  the 
Select  Cellstring  tool  and  then  double  click  on  the  box  that  lies  in  the  center  of  the  new  cell  string. 
This  action  will  activate  the  M2D  Boundary  Conditions  dialog  shown  in  Figure  2.  Click  on  WSE- 
forcing  on  the  left  side  of  the  dialog.  This  action  enables  the  user  to  choose  the  type  of  water-level 
forcing  information  that  M2D  will  apply.  Click  on  Extract  from  data  set  and  then  click  the  Extract 
BC  button  to  bring  up  the  Extract  Boundary  Conditions  dialog.  There  are  two  areas  of  this  dialog, 
with  the  top  portion  being  the  Time  Control  and  the  bottom  portion  being  the  Select  Dataset  dialog. 
The  Time  Control  lets  the  user  select  a  portion  of  the  available  solution  that  can  be  extracted.  The 
Select  Dataset  provides  a  display  of  the  data  sets. 

For  Time  Control,  the  SMS  will  automatically  set  the  information  extraction  to  take  place  from  the 
beginning  of  the  available  solution  file.  The  length  of  the  information  extracted  (in  time)  will 
correspond  to  the  specified  duration  of  the  M2D  simulation.  The  user  can  select  other  portions  of 
the  solution  if  needed  by  clicking  on  Regional  Model,  then  editing  the  starting  and  ending  times. 
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To  extract  water  levels  from  an  ADCIRC  solution,  select  water-surface  elevation  (63)  as  shown  in 
Figure  3  (the  (63)  refers  to  the  default  file  extension  for  global  water-surface  elevation  files  written 
by  ADCIRC).  This  action  will  specify  the  solution  from  which  water-level  values  will  be  obtained. 
Click  Extract  to  perform  the  extraction  from  the  ADCIRC  solution  and  interpolation  to  the  cell 
centers  of  the  boundary  cells  for  M2D. 

After  the  extraction  is  complete,  the  user  can  view  the  water-level  curves  by  clicking  View  in  the 
M2D  Boundary  Conditions  dialog.  The  Extracted  Data  dialog  will  open  and  a  blank  plot  will  be 
shown  along  with  a  list  of  cells  for  which  data  have  been  extracted  and  interpolated  (Figure  8).  Time 
series  of  water  levels  can  be  viewed  for  selected  boundary  cells  by  clicking  in  the  boxes  under  View, 
as  shown  in  Figure  9.  Plot  options  can  be  modified  by  clicking  on  the  plot,  then  right  clicking. 
Click  OK  to  exit. 


Figure  8.  Initial  Extracted  Data  dialog 

For  this  example,  water  levels  from  ADCIRC  have  also  been  extracted  at  the  Shinnecock  Canal, 
which  is  located  on  the  northernmost  section  of  Shinnecock  Bay.  This  canal  has  a  tidal  gate, 
allowing  water  to  enter  Shinnecock  Bay  but  not  leave  it. 

After  the  extracted  boundary  conditions  and  all  other  M2D  settings  have  been  specified,  the  project 
can  be  saved  by  clicking  on  File  at  the  upper  left  portion  of  the  screen,  then  Save  Project.  M2D  can 
then  be  launched  by  clicking  on  the  M2D  menu  at  the  top  of  the  screen,  then  Run  M2D  from  the 
drop-down  menu.  If  the  SMS  cannot  find  the  M2D  executable,  the  user  will  need  to  specify  where 
the  model  is  located  on  the  computer. 
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Figure  9.  Extracted  Data  dialog  with  two  water-level  curves 

Selected  results  of  the  M2D  simulation  at  ebb  and  flood  tide  are  provided  as  an  example  solution. 
Water  level  contours  and  current  velocity  vectors  are  shown  in  Figure  10  during  ebb  tide.  White 
areas  in  the  bay  denote  dry  cells.  This  snapshot  was  taken  near  the  end  of  an  ebb  cycle  and  the  tidal 
current  in  the  coastal  ocean  has  started  to  reverse. 

Current  speed  contours  and  direction  vectors  during  flood  at  the  inlet  and  bay  are  shown  in 
Figure  11.  The  dry  areas  of  the  bay  are  smaller  than  those  in  Figure  10,  indicating  higher  water 
level.  Strong  currents  are  formed  within  the  inlet  and  at  the  Shinnecock  Canal.  In  this  snapshot,  the 
coastal  tidal  current  is  flowing  from  east  to  west.  This  flow  pattern  along  the  coast  can  only  be 
obtained  by  specification  of  varying  water  levels  along  the  ocean  boundary. 

DISCUSSION :  This  CHETN  has  described  coupling  of  the  regional  and  local  circulation  models 
ADCIRC  and  M2D  within  the  SMS.  This  type  of  coupling  is  advantageous  for  local  applications  in 
which  multiple  simulations  are  required.  The  SMS  automates  the  model  coupling  by  extracting 
water-level  values  from  global  ADCIRC  solutions,  interpolating  them  to  M2D  boundary  cell  centers, 
and  saving  the  time  series  as  M2D  forcing  input  fdes.  If  solutions  from  previous  ADCIRC 
simulations  are  available,  the  user  may  not  have  to  run  the  regional  model,  but  can  apply  existing 
solutions  as  boundary  conditions. 
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Figure  10.  Water  level  and  current  vectors  at  ebb  tide  over  entire  domain 


Figure  1 1 .  Current  speed  and  direction  at  flood  tide  over  inlet  and  bay 
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Coupling  of  ADCIRC  and  M2D  has  been  conducted  for  Shinnecock  Inlet,  NY,  and  at  Bay  Center, 
WA  (Militello  2002;  Militello,  Parry,  and  Arden  2002).  Shinnecock  has  a  tide  range  of  about  1.5  m 
(5  ft)  and  a  flood  shoal  that  becomes  inundated  and  exposed  during  the  tidal  cycle.  Bay  Center  is 
located  in  Willapa  Bay,  which  has  a  tide  range  of  about  2.7  m  (9  ft).  Approximately  half  of  Willapa 
Bay  dries  during  low  tide.  In  both  cases,  the  water-level  values  transferred  from  ADCIRC  to  the 
local  M2D  model  specified  high-quality  boundary  conditions  and  included  tidal  phase.  This 
boundary  specification  enabled  accurate  representation  of  processes  at  both  locations. 

ADDITIONAL  INFORMATION:  Questions  about  this  CHETN  can  be  addressed  to  Dr.  Adele 
Militello  at  (707)  444-2173  or  CoastalAnalysis@cox.net ,  to  Dr.  Alan  K.  Zundel  at  (801)  378-2812  or 
Zundel@byu.edu ,  or  to  Ms.  Mary  A.  Cialone  at  (601)  634-2179  or  Mary.A.Cialone@erdc. 
usace.army.mil.  This  CHETN  was  produced  under  the  Inlet  Modeling  System  work  unit  of  the 
Coastal  Inlets  Research  Program.  For  information  about  the  Coastal  Inlets  Research  Program,  please 
contact  the  Program  Manager,  Dr.  Nicholas  C.  Kraus,  at  (601)  634-2016  ox  Nicholas.  C.Kraus@erdc. 
usace.army.mil.  This  technical  note  should  be  cited  as  follows: 

Militello,  A.,  and  Zundel,  A.  K.  (2002).  “Coupling  of  regional  and  local  circulation 
models  ADCIRC  and  M2D,”  ERDC/CHL  CHETN-IV-42,  U.S.  Army  Engineer 
Research  and  Development  Center,  Vicksburg,  MS. 
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